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Vertically pointing Doppler radar has been used to study 
the evolution of ice particles as they sediment through a cir- 
rus cloud. The measured Doppler fall speeds, together with 
radar-derived estimates for the altitude of cloud top, are 
used to estimate a characteristic fall time t c for the 'average' 
ice particle. The change in radar reflectivity Z is studied as 
a function of t c , and is found to increase exponentially with 
fall time. We use the idea of dynamically scaling particle size 
distributions to show that this behaviour implies exponen- 
tial growth of the average particle size, and argue that this 
exponential growth is a signature of ice crystal aggregation. 



1. Introduction 

The growth of ice crystals and aggregate snowfiakes in 
clouds is a key process both for the development of precip- 
itation (Jiusto and Weickmann 1973), and in terms of the 
effect such clouds have on climate (Houghton 2001). In this 
work, we use radar observations of deep cirrus to study the 
growth of ice particles as they sediment through the cloud. 

Vertically-pointing measurements of radar reflectivity Z 
and Doppler velocity Vd were made using the 35 GHz 
(8.6 mm) 'Copernicus' radar at the Chilbolton Observatory 
in southern England. At this wavelength, the overwhelm- 
ing majority of cirrus-sized particles are within the Rayleigh 
regime where the backscattered intensity is proportional to 
the square of particle mass m: 



longer a detectable radar return, we calculate the fall time 
associated with height h as: 



Z = 



36 \K ic 



0.93tt 2 p 



x / n(m)m 2 dm, 



(1) 



where p lcc is the density of solid ice and n(m)dm is the num- 
ber density of particles with mass between m and m + dm. 
The dielectric factor Ki cc contains the information about the 
shape and dielectric strength of the particles: for spherical 
ice particles Ki ce = and the permittivity of ice e at mil- 
limetre wavelengths is approximately 3.15 (Jiang and Wu 
2004). The Rayleigh scattering approximation at 35 GHz 
is accurate to within 10% for particles with a maximum di- 
mension of 1 mm or less (Westbrook et al 2006). 

The Doppler velocity is Vd — v t + fair, where v t is the 
m 2 -weighted average terminal velocity of the particles and 
Vait is the vertical air motion. We use these measurements 
to estimate a characteristic particle fall time t c , which we 
define as the time for which the 'average' particle (with ter- 
minal velocity Vt) has been falling. Note that the Doppler 
velocity is weighted by the reflectivity making it sensitive 
to the larger ice particles, and so our average fall time will 
also be weighted toward these large particles. Taking the 
cloud top height h t0 p to be the altitude at which there is no 
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Given this new measure, we are in a position to investigate 
the evolution of the ice particles, by studying the variation 
of reflectivity Z with increasing fall time t c . The advan- 
tage of this method, as opposed to simply studying Z as 
a function of height, is that t c represents the physical time 
for which the average ice particle has been falling to reach 
a given height h, allowing us to relate our results to the- 
oretical models of ice particle growth. Note that we have 
implicitly assumed that the cloud is in a steady state, such 
that the properties of the ice particles at height h do not 
change significantly over the length of time it takes a parti- 
cle to fall from cloud top to cloud base (which is between 45 
minutes and 2 hours for the cases shown here). Essentially 
this means that the cloud does not evolve significantly on 
this time scale and is advecting as a rigid body across the 
radar beam. We therefore apply our technique only to non- 
precipitating, well developed ice clouds where there is there 
is low wind shear. 

2. Cloud Data 

Our case study is a cirrus cloud observed over Chilbolton 
on the 13 th of May 2004. The temperature at cloud top (as 
forecast by the Met Office mesoscale model, Cullen 1993) 
was approximately — 40° C, and the cloud base was close to 
— 15°C; the average wind shear over the depth of the cloud 
was approximately 2ms _1 km _1 . Measurements of reflec- 
tivity and Doppler velocity were made and the time series 
of these observations is shown in figure 1. The radar gate 
length is 30 m (Illingworth et al 2006). The values of Z 
and Vd are averages over periods of 30 seconds: in figure 
lc we also show the standard deviation a— of the 1-s aver- 
age Doppler velocity over each 30-s period, to indicate the 
small-scale variability in Vd- This measure allows the level of 
turbulence in the cloud to be assessed (Bouniol et al 2003). 

Figure 2 shows four representative vertical profiles sam- 
pled from different portions of the cloud, indicated by the 
dashed lines on figure 1. Ten consecutive 30-s profiles were 
averaged over a period of ~7 minutes in order to smooth out 
the variability caused by fall streaks in the data. The high- 
est detectable cloud pixel (corresponding to ~ —15 dBZ) 
from the profile is taken as a measure of cloud top. The fall 
time at each height bin is calculated from the Doppler ve- 
locity profile as per equation 2, and we plot Z as a function 
of t c . From figure 2 we see that reflectivity increases rapidly 
with fall time (note the logarithmic dBZ units), which we 
interpret as rapid growth of the ice particles. This could 
potentially be occurring through a number of possible mech- 
anisms: deposition of water vapour; aggregation via differ- 
ential sedimentation of the ice particles; or collisions with 
supercooled drops (riming). In section 4 we show that it 
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is likely that aggregation is the dominant growth mecha- 
nism. The increase in Z appears to be exponential to a 
good approximation, and occurs for between 2500 and 5000 
seconds in the profiles shown here. The slopes on the log 
scale vary between approximately 2.5 x 10 dBZ s _1 and 
5 x 10~ 3 dBZ s _1 , presumably depending on how much ice 
is being produced at cloud top. After this time there is a 
sharp turn over in the Z(t c ) curves, and we attribute this to 
evaporation of the particles near cloud base. Such evapora- 
tion often results in increased air turbulence for which the 
particles themselves act as tracers, resulting in large vari- 
ability in the Doppler velocity. In the earlier profiles (07:09 
and 07:38 UTC) this was not evident; however, in the later 
profiles (08:06 and 08:30 UTC) the higher ice water content 
and time-integrated evaporative cooling triggered convective 
overturning and turbulence, and this is reflected in our ob- 
servations (figures 1 and 2), which show a sudden increase 
in a— at approximately the same time as the turn over in 
Z(t e ). 

Exponential growth has also been observed in a number 
of other cloud data sets, and four more example profiles from 
well developed non-precipitating ice clouds during April and 
May 2004 are shown in figure 3. This is an interesting fea- 
ture of the data, and a robust one in the face of errors in 
ht op : if Z(t c ) is exponential, then even if we have underes- 
timated the cloud top somewhat (on account of the limited 
sensitivity of the radar), this will merely correspond to an 
offset in the fall time, and the exponential shape of Z(t c ) is 
still preserved. It is interesting to note that the transition 
from growth to evaporation is not always sharp as it is for 
the 13 th May profiles: we speculate that this may be the 
result of aggregation continuing to some extent within the 
evaporation layer. 

3. Scaling analysis 

Here we show how the reflectivity Z is related to the av- 
erage particle size. Scaling or 'normalised' forms for the size 
distributions of both liquid and ice particles have been pro- 
posed in a number of recent articles (rain: Testud et al 2001, 
Illingworth and Blackman 2002, Lee et al 2005; ice: Field 
and Heymsfield 2003, Westbrook et al 2004a, b, Delanoe et 
al 2005). The essence of these rescaling schemes is that the 
underlying shape of the distribution cj>(m/{m}) is the same 
throughout the vertical profile, but is rescaled as a function 
of the (increasing) average particle mass (m) as the particles 
grow: 



n(m) = IWC x {my 2 4>{m/{m)) . 



(3) 



where we have normalised by the ice water content IWC. 
The universal function <f> is dimensionless. Equation 3 in- 
dicates that a single average particle mass (m) is sufficient 
to characterise the evolution of the particle size distribution 
(relative to the IWC or some other moment of the distribu- 
tion), and this is key to our analysis. 

An example of such a distribution is that assumed in the 
UK Met Office's Unified Model (Wilson and Ballard 1999). 
Mass m and diameter D are assumed to be in a power law 
relationship m = a' D b , with an exponential distribution for 
particle diameter: 



n'(D) = N exp(-AD), 



(4) 



where n'(D) — n(m)dm/dD. A single bulk prognostic vari- 
able is used for the ice particle mixing ratio and No is param- 
eterised to decrease with increasing temperature to mimic 
particle growth. The parameter A is calculated from the 
predicted IWC and No, and is interpreted as a reciprocal 



average diameter (D) 1 (eg. Brown et al 1995). Within the 
framework (3) above, this distribution corresponds to: 



t/>{x) = [bT(b + l)] 



-i (i- 



b)/b 



exp (-x 1/b ) , 



(5) 



where x = m/(m), and (m) = a'(D) b . 

Irrespective of what form is assumed for (f>(x), a scaling 
relationship between different moments of the distribution 
may be found. The k th moment of the mass distribution is 
given by: 



poo poo 

/ n{m)m k dm = (m) fc_1 IWC / </>(x)x k dx. (6) 
Jo Jo 



Note that <j)(x)x k dx is a dimensionless constant. Simi- 
larly, the radar reflectivity (1) is given by: 



Z = (m) IWC 



36\K icc [ 
0.9371-^2 
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(x)x dx. 



(7) 



Combining these two equations we may relate Z to an arbi- 
trary moment Mk of the distribution: 



Z = (mf- k [ M k x 
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At this point we make a crucial assumption: that there is 
some moment of the distribution k which is approximately 
constant through the vertical profile. In the case where ag- 
gregation is the dominant growth mechanism with a fixed 
production of ice mass at cloud top, one would expect the 
mass flux density of ice n(m)mv(m)drn to be constant. 
Mitchell (1996) indicated that a power law for ice particle 
fall speeds v is a good approximation: v(m) oc m c , so for 
pure aggregation k = 1 + c. Similarly, where diffusional 
growth or riming is dominant, the total number flux of par- 
ticles would be roughly constant and k — c would be the 
conserved moment. If this assumption holds then the brack- 
etted expression (...) in equation 8 is fixed through the ver- 
tical profile, and Z oc (m) 2 ~ fe . Given our observations of 
exponential Z(t c ) and the predicted power law between Z 
and (m) above, we conclude that the average particle mass 
is growing exponentially with fall time. 

4. A Signature of Aggregation? 

We offer a possible explaination for the exponential 
growth of ice particles described above. Aircraft observa- 
tions have indicated that aggregation is often the dominant 
growth mechanism for particles larger than a few hundred 
microns in cirrus clouds (Field and Heymsfield 2003) , and it 
is these large particles which dominate the radar reflectivity. 

Recently, Westbrook et al (2004a,b) modelled ice particle 
aggregation by considering a rate of close approach between 
pairs of ice particles with masses m and m! : 



K = j (D max + D' max ) 2 \v - v'\ , 



(9) 



where v and D max are the associated fall speed and maxi- 
mum dimension. Particles were picked according to the rate 
above, and traced along possible trajectories to accurately 
sample the collision geometries of the non-spherical ice par- 
ticles. The fall speeds were prescribed in the vein of Mitchell 
(1996): 



(10) 
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where the adjustable parameter a determines the drag 
regime (inertial flow a = i; viscous flow a = 1). One of the 
key results from these simulations was that the aggregates 
produced by the model had a power law relationship between 
mass and maximum dimension m oc D b max , where the expo- 
nent is determined purely by the drag regime: b = 1/(1 — a) 
for a < |. This relation is also backed up by a theoret- 
ical argument based on a feedback between the aggregate 
geometry and collision rate (Westbrook et al 2004b). For 
large snowflakes, a — > ^ and b — > 2, in good agreement with 
aircraft observations (eg. b = 1.9, Brown and Francis 1995; 
b = 2.04, Heymsfield et al 2002). 

In this study we are interested in the average ice particle 
growth rate, which is determined through the scaling of the 
collision kernel (9). Given the above relationship between a 
and b, and equations 9 and 10, we see that if one doubles 
the masses of the aggregating particles m,m': 

K(2m,2m') = 2 x K(m,m), (11) 

where A = a + 1/6 = 1. This parameter controls the scaling 
of the particle growth rates and as such controls the growth 
of the average particle mass. 

Van Dongen and Ernst (1985) have shown that the co- 
agulation equation (Pruppacher and Klett 1997) has solu- 
tions with the same scaling form as (3), and predicts that 
the average particle mass grows according to the differential 
equation: 



where w is a constant. Given our prediction of A = 1 from 
the aggregation model: 

(m) oc exp(wt c ), (13) 

i.e. the prediction from aggregation theory is that average 
particle mass grows exponentially with fall time, in agree- 
ment with our observations. We note that the Van Don- 
gen and Ernst analysis is for cases where total mass is con- 
served: however given the observed scaling behaviour (3) 
and a power law relationship between mass and fall speed, 
the case where mass flux density is conserved should yield 
the same result. 

The growth of particles by diffusion of water vapour may 
also be described by a similar equation to (12). However 
in that case A = 1/6 and w = 4nCoa~ 1/b x s/(A + B), 
where Co is the 'capacitance' per unit diameter, s is the 
supersaturation with respect to ice, and the terms A and 
B depend on temperature T and pressure P (Pruppacher 
and Klett 1997). For a given set of conditions (s,T, P), the 
growth by deposition would be expected to increase slower 
with particle size than for aggregation, taking a power law 
form (m) oc t b J < ^ b ~ V> . In real clouds these conditions do not 
stay constant, and there is a correlation between increasing 
particle size and increased temperature and supersaturation, 
which could lead to a faster growth rate. However, it would 
take a considerable conspiracy between these variables to 
obtain a constant exponential growth throughout such an 
extensive region of the cloud as is observed in our radar 
data. It also seems extremely unlikely that this correlation 
would be the same for all five cirrus cases shown in figures 
2 and 3. We note that there is a region of sub-exponential 
growth close to cloud top (small t c ) in some of the profiles 
in figure 3: we suggest that it is in this region, where the 
particles are small and falling slowly, that diffusional growth 
dominates. 



It seems very unlikely that riming dominated the ice par- 
ticle growth: a large number of supercooled drops through- 
out the depth of the cloud would be required for this to be 
the case. Given the cold temperatures in the cloud (between 
— 40°C and — 15°C as discussed in section 2), it is very un- 
likely that supercooled drops would persist on long enough 
time scales and in large enough quantities to dominate the 
growth over the 2.5 km or so for which we have observed 
Z(t c ) to increase exponentially. We therefore discount de- 
position and riming, and assert that our observations are an 
indicator that aggregation is the dominant growth mecha- 
nism for the ice particles in these clouds. 

5. Discussion 

Doppler radar measurements of cirrus cloud were used to 
study the evolution of the ice particles sedimenting through 
it. The results indicate that in the cases studied the av- 
erage ice particle mass grows exponentially with fall time, 
in agreement with the theoretical expectation for aggrega- 
tion, and we believe that this is evidence that aggregation 
of ice crystals is the dominant growth mechanism for large 
particles in deep, well developed ice clouds. 

Vertical profiles of reflectivity in ice have been much stud- 
ied in order to estimate rainrates at the ground. Fabry 
and Zawadzki (1995) observed an approximately constant 
d(dBZ)/d/i gradient, and used this to rule out deposition 
as a growth mechanism. This may be linked to our cirrus 
observations; however their results were near the melting 
layer, and Z was ~20dB higher. We have compared profiles 
of dBZ-/i and dBZ-£ c for our cirrus cases and find that while 
the dBZ-£ c profiles are straight lines with a constant gradi- 
ent, the dBZ-/i profiles have an appreciable curve to them. 
The fact that our analysis 'straightens' these curved profiles 
is good evidence that our approach of using the Doppler 
velocity to estimate t c from h is an appropriate one, and 
that aggregation is controlling the distribution of large ice 
particles. 

The constant w described in the aggregation theory above 
is directly related to the mass flux density, so measurements 
of the dBZ-t c slope may allow the derivation of this quantity, 
and the data could be combined with Doppler velocity mea- 
surements to estimate the ice water content. However, the 
sticking efficency of the ice particles (which we assume to be 
constant with particle size) is also a factor in to, and this is a 
parameter for which there are few reliable experimental es- 
timates. For warmer, 'stickier' ice crystals at temperatures 
above — 5°C this may be more feasible since the sticking 
efficiency should be close to unity. 

We have assumed the ice particles fall vertically. In re- 
ality there is likely to be some horizontal shear, and this, 
combined with variability in ice production of the cloud-top 
generating cells results in visible fall streaks (see fig. 1). 
Size-sorting along the streaks (Bader et al 1987) is a poten- 
tial source of error in our analysis; however, by averaging 
the reflectivity profiles over ~7 minutes of data we have 
been able to ameliorate it considerably. 

Directions for future work are to make dual-wavelength 
radar measurements of cirrus in order to obtain a more di- 
rect estimate of particle size (Westbrook et al 2006). This 
would help to pin down the dominant growth mechanism, 
allowing us to study moments other than Z, and analyse 
whether k — 1 + c (aggregation) or k = c (deposition, rim- 
ing) is the moment conserved through the cloud. Aircraft 
observations (Field et al 2005) have indicated a broadly ex- 
ponential trend between Z and temperature - it would be 
valuable to combine simultaneous radar and aircraft mea- 
surements to see if the exponential growth in Z with t c is 
accompanied by exponential growth in (m) and increased 
concentrations of aggregates. Also, further studies of other 
cirrus cases, both at Chilbolton and other radar sites, could 
be of interest to see how widespread the observed exponen- 
tial trend is. 
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Figure 1. Radar time series of a cirrus cloud over 
Chilbolton on the 13 th May 2004. Panels (a) and (b) 
show the reflectivity Z and Doppler velocity Vd respec- 
tively: both are averages over 30 s of data. Panel 
(c) shows the standard deviation a- of the 1-s average 
Doppler velocity for each 30-s period, indicating the vari- 
ability in v d . 
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Figure 2. Four 'snapshot' vertical profiles from the cir- 
rus case, taken at 07:09, 07:38, 08:06, and 08:30 UTC. 
Each profile shown is the average of ten consecutive 30-s 
profiles. Top row is reflectivity in dBZ as a function of 
characteristic fall time (points) . The solid line is intended 
to guide the eye, and indicates an exponential growth in 
Z with t c . Bottom row is a— as a function of t c , which 
we use as an indicator of particle evaporation near cloud 
base. 
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Figure 3. Vertical profiles of reflectivity from a sample 
of four more cirrus cases measured over Chilbolton during 
April and May 2004. Exact times and dates are indicated 
on the individual panels. All show an exponential growth 
of reflectivity with fall time over a significant portion of 
the cloud vertical profile. 



